Analysis of the exciton-exciton interaction in semiconductor quantum wells 
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The exciton-exciton interaction is investigated for quasi- two-dimensional quantum structures. 
A bosonization scheme is applied including the full spin structure. For generating the effective 
interaction potentials, the Hartree-Fock and Heitler-London approaches are improved by a full two- 
exciton calculation which includes the van der Waals effect. With these potentials the biexciton 
formation in bilayer systems is investigated. For coupled quantum wells the two-body scattering 
matrix is calculated and employed to give a modified relation between exciton density and blueshift. 
Such a relation is of central importance for gauging exciton densities in experiments which pave the 
way toward Bose-Einstein condensation of excitons. 
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I. INTRODUCTION 

Excitons (bound pairs of electron and hole) play a 
key role in semiconductor optics, and their theoretical 
understanding has a long history. 1 Although excitons 
obey Bose statistics, the fermionic structure of their con- 
stituents is always important and forbids to treat excitons 
as nearly ideal Bose particles. Among the first attempts 
for bosonization of excitons, we quote Refs. 2-4. Both 
the fermionic exchange and the Coulomb forces between 
two excitons can be condensed into an effective exciton- 
exciton (XX) interaction potential. Efforts for deriv- 
ing XX potentials are abound in the literature, cover- 
ing three-dimensional excitons in bulk semiconductors as 
well as quasi- two-dimensional excitons in quantum wells. 
However, almost all of these attempts 5-8 were restricted 
to the Hartree-Fock level, i.e., taking into account the XX 
interactions up to second order in the elementary charge, 
e 2 . Sometimes it has been claimed that this would be 
enough for treating sufficiently accurate the exciton gas 
at low density nx- This is not correct since already at 
zero density, the single-exciton bound state calls for an 
infinite summation in powers of the Coulomb potential, 
and consequently the next term (linear in nx) cannot be 
truncated either. It was not before 2001 that in a semi- 
nal paper by Okumura and Ogawa 9 the first XX poten- 
tial beyond Hartree-Fock 10 for bulk semiconductors has 
been constructed, in close analogy to the Heitler-London 
approximation in atomic physics. 

The interest in a proper description of XX interac- 
tions has been surely intensified by the actual search for 
Bose-Einstein condensation of excitons, which has been 
predicted theoretically already decades ago. 11-13 Due to 
the small exciton mass compared to atomic systems, the 
critical temperature for the condensate is expected to be 
a few kelvins for a density of 10 17 cm -3 in bulk GaAs, 
within easy reach experimentally A fundamental prob- 
lem, however, is the finite life time of the excitons, which 
hinders the relaxation into thermal equilibrium. One pos- 
sible way out are coupled quantum wells (CQWs) which 
came into focus a few years ago. 14-18 A static electric 



field in growth direction forces electrons and holes to re- 
side in adjacent quantum wells which are separated by a 
barrier. Due to this spatial separation, the indirect exci- 
tons exhibit extremely long life times, which is a good 
condition for reaching thermal equilibrium. However, 
these spatially indirect excitons form dipolcs leading to 
a strong and long-range repulsion, which complicates the 
theoretical description as well as the experimental real- 
ization of a dense cold exciton gas. 19 ' 20 Recently large 
progress has been reported for spatially indirect excitons 
in electrostatic 21 as well as optical traps. 22,23 

These practical demands led us to investigate the XX 
interaction in more detail, with special emphasis on cou- 
pled quantum wells. In particular, we improve on the 
Heitler-London-type treatment in Ref. 9 by solving the 
four-particle Schrodinger equation for two electrons and 
two holes numerically. Before doing so we address the 
complex spin structure of the exciton composed of a spin 
1/2 electron and a spin 3/2 heavy hole (Sec. II). The im- 
portance of spin-dependent effects in the exciton gas has 
been emphasized e.g., in Ref. 24. In the limit of immo- 
bile holes (their mass being usually much larger than that 
of the electrons), we can derive effective spin-dependent 
XX potentials (Sec. Ill) which contain in addition to the 
dipole-dipole repulsion and exchange effects the weak van 
der Waals forces. The latter is of importance to rectify a 
recent claim 25 ' 26 to have explained the bead pattern for- 
mation which appears at low temperatures in the ring- 
shaped CQW emission. 27-29 

With a proper XX interaction at hand, we are able 
to calculate biexciton states (excitonic molecules), which 
are well known from bulk semiconductors and have 
been observed in single quantum wells. 30,31 Since biex- 
citon energies have been calculated with high precision 
elsewhere, 32 we use these results to judge our approxi- 
mate treatment in Sec. IV. For the CQW situation sim- 
plified to a bilayer system, we identify the parameter val- 
ues (essentially mass ratio and charge separation) which 
limit the existence of biexcitons. 

Addressing the many-exciton case, we will treat the 
excitons as effective bosons with a renormalized interac- 
tion potential, derived from the underlying electron-hole 



description. 33 As an application, we investigate in Sec. V 
two realistic CQW structures. The numerically gener- 
ated XX potentials are used to calculate two-exciton scat- 
tering phase shifts which are the main ingredient for a 
T-matrix based quasi-particle dispersion 34 . For the low- 
density case, we are able to calculate excitonic blueshift 
and scattering- induced broadening linear in nx (Sec. VI) . 
We find a stunning reduction of the blueshift compared 
to the simple "capacitor formula" and relate this finding 
to features in the XX pair-correlation function which - 
due to the strong repulsion - resembles more a Fermi gas 
than a free Bose gas (Sec. VII). 

Conclusions are drawn in Sec. VIII, while a few tech- 
nical details are deferred to the Appendix. 

II. MANY-EXCITON HAMILTONIAN 

To derive the many-exciton Hamiltonian we follow the 
work of de-Leon and Laikhtman. 8 ' 35 The task is to find 
the matrix elements of the electron-hole Hamiltonian 
with an appropriate two-exciton wave function and to 
implement them into an effective bosonic Hamiltonian 
for many excitons. The effective-mass Hamiltonian for 



two electrons and holes reads 

h 2 A el h 2 A hl h 2 A e2 h 2 A h2 

2ch 2m e 2m/, 2m e 2m/, 
+ v ee (r el - r e2 ) + v hh (r hl - r h2 ) - v eh (r el - r hl ) 
- v eh (r e2 - r h2 ) - v eh (r el - r h2 ) - v eh (r e2 - r M ), (1) 

where the interaction part [second and third line of 
Eq. (1)] is composed of Coulomb interactions between 
particle a and b, which can be either electron (a = e) 
or hole (a — h), 

e 2 e 2 

v ab (r) = -± with e 2 = - . (2) 

r 47re e s 

Please note that the corresponding sign (attractive or 
repulsive) has been made explicit in Eq. (1). 

We define exciton quantum field operators ^|(R) 
which create an exciton at exciton center-of-mass (cm.) 
position R with spin s, and \I> a (R), which annihilates the 
same exciton. The subsequent Hamilton operator will be 
written as 



h xx = f d R^M/t( R) _!^M, s (R) + I j d ndn' ^w,.3rf(R-R')*Ii(R)*l2(B-')*d(R')*riW. (3) 

J ■ J sls2s3s4 



with the exciton mass M = m e + m/, in the kinetic 
energy. W s i S 2,«3a4(R) is the spin and space-dependent 
pair interaction potential. In the following, it will be 
extracted from a careful study of the four-particle prob- 
lem [Eq. (1)]. The exciton spin index s = s e + Jh is the 
sum of the electron spin (s e = ±1/2) and the heavy- 
hole angular momentum (Jh = ±3/2). We neglect the 
light holes, which are separated due to the confinement 
effects in the quasi- two-dimensional quantum well (QW). 
Now, s runs over four values s = ±1 (bright states) 
and s = ±2 (dark states). Since electron and hole are 
fcrmions, we have to use a properly symmetrized ansatz 
for both component's spin wave functions y^ e (s e \,s e2 ) 
and xl/,(4ii42), which make the overall wave- func- 
tion antisymmetric with respect to the exchange of the 
electrons and the holes, respectively. The labels p and 
q denote the parity which can be symmetric (p = s) or 
antisymmetric (p = a). Together with the spatial part 
we write the total two-exciton wave function, 

*£(r e ii Sel,r e2 , S e2 , r hu J h i,r h2 , J h2 ) 

= ^(rei,r e2 ,r hl ,r h2 )Xee P (sei,s e2 )xhh(Jhi, Jh2), (4) 

where the upper label stands for the electron part, and 
the lower one for the hole part. Please note that due to 



the fermionic nature of the particles, the parities of the 
spin part and of the spatial part have to be opposite to 
each other. Taking now the matrix elements of Eq. (1) 
with the ground state wave functions we get four different 
potentials, 36 

UZ = {W q \H% h \W q ). (5) 

Due to the one to one correspondence between the ex- 
citon spin and the spin of its constituents, it is possible 
to express the basis vectors in the space of symmetric 
and antisymmetric spin wave functions by the spin eigen- 
states of the excitons |si,S2). After a straightforward 
unitary transformation one gets the interaction matrix 
elements W s i S 2. S 3s4(R) as shown in Table I for columns 
| S3S4) and rows {s x s 2 \. Here Uf = 1/2 (C/ g a ± U s q ) and 
similar for the lower hole index. We can clearly see 
the block structure of the interaction which reflects the 
conservation of the total spin si + s 2 = s 3 + s 4 . The inter- 
action channels can be classified as follows: There is the 
direct channel where an initial state | si, s 2 ) will remain 
unchanged. The other channel is of exchange type with a 
change in the initial state due to three different processes: 
Exchange of electrons, exchange of holes, and exchange 
of both simultaneously. The electron-hole exchange pro- 
cess (longitudinal-transverse splitting of the exciton) has 



2 





2,2) 


|1,1> 


-1,-1) 


1-2,-2) 


|2,-1) 


1-1,2) 


|l,-2) 


-2,1) 


{ 2, 2 1 


U° 
















(Ml 




u<< 










n 
u 




















(-2, -2 








Ui 










(2,-11 
(-1,2 











yo 


U" 






(1, -2| 












Ui 




(-2, 1 














ut 


!/" 




1,2} 


|2,1) 


-2,-1) 


1-1. -2) 


|2,-2) 


1,-1) 


1-2,2) 


1-1,1) 


< 1, 2 1 


ui 


v- 














(2,1 


ij- 


ui 










n 
u 




(-2, -1 
















(-1, -2 






u~ 


ui 










(2, -2 














u- 




(1,-1 
(-2,2 
(-1,1 













u' 

UZ 


u+ 

ui 

ui 


u- 
Vi 

ui 


Table I 


Spin structure of the exciton-exciton interaction 


WsU2,s3s4,(R) 


in quantum wells. 
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Table II: Spin structure of the exciton-exciton interaction 
W s is2,s3s4(R) in the case of infinitely heavy holes. 



been neglected already in the starting Hamiltonian. 

In a case where the hole mass is much larger than the 
electron mass, the exchange of holes becomes negligible 
and the four different potentials [Eq. (5)] collapse into 



two, U a 



U° and U s 



Ug, which correspond to wave 



functions properly symmetrized only with respect to the 
electrons. In this case the interaction part of the Hamil- 
tonian simplifies to the form shown in Table II where 
U = l/2(U a ±U S ) and can be cast into the form given 
in the Appendix, [Eq. (Al)]. 



III. EFFECTIVE INTERACTION POTENTIALS 

To derive the spatial dependence of the interaction po- 
tentials introduced in Sec. II, approximations have been 
introduced in the literature. In the infinitely heavy- 
hole limit, the Hcitlcr-London ansatz is well known from 
atomic physics and has recently been brought into exci- 
ton physics by Okumura and Ogawa. 9 In a bulk semicon- 
ductor, the two excitons resemble a hydrogen molecule in 
this limit. The problem of four particles (two electrons 



and two holes) simplifies here to a two-particle problem 
for the electrons, while the hole-hole distance R enters 
as a parameter. Therefore, the Coulomb potential be- 
tween the holes gives just a fixed additional term in the 
Hamiltonian. The two-exciton wave function can then be 
written as a properly antisymmetrized product of single- 
exciton wave functions in the Is ground state </>(r), cen- 
tered around the position of each hole. For the spatial 
part of the two-exciton wave function, we can write 



^ s ' a (ri,r 2 ) 



1 0(n)0(r a -R) ± 0(n-R)^(r 2 ) 



V2 V 1 ± ^ 2 ( R ) 

with the wave function overlap, 



O(R) = J dr0(r)0(r-R), 



(G) 



(7) 



which leads to the potentials, 

U d (R) - U X (R) 



II s (R) 
U a (R) 



1 + C 2 (R) ' 
U d {R) + U X (R) 
1-0 2 (R) ' 



(8) 



with direct Ud(R) and exchange U X (R) potential [ex- 
plicit expressions are given in the Appendix, Eq. (A4) 
andEq. (A5)]. 

A simplified version of the Heitler-London approx- 
imation, which has usually been applied to excitonic 
systems, 5,6,8,37 is the famous Hartree-Fock treatment 
where the normalization denominator in Eq. (6) is left 
out and is consequently missing in Eq. (8). This scheme 
can easily be generalized to arbitrary hole masses. How- 
ever, it leads to a nonlocal exchange potential 5 and to 
problems with the orthogonality of the basis states. 38 

To improve over the approximations discussed so far, 
we have solved the Schrodinger equation for the two elec- 
trons numerically using the Lanczos algorithm, treating 
the holes as infinitely heavy and thus immobile. The 
equation to be solved reads 



— (A ri + A r2 ) + v ee (r 2 - ri) + v hh (R) - v e h(ri) 
2/i 



-v eh (r 2 ) - v eh (r 2 - R) - v eh (r 1 - R) 



* p (ri,r 2 ) 



= (-2B x + U p (R))W(r 1 , 



ra), 
(9) 



where fj, is the effective electron mass in the plane and 
Bx is the single-exciton binding energy. The zero of en- 
ergy is chosen to be the band gap. Fig. 1 shows the re- 
sulting effective interaction potentials in Hartree-Fock, 
Heitler-London, and full numerical quality exemplarily 
for a strictly two-dimensional system. We note the un- 
physical behavior of the Hartree-Fock antisymmetric po- 
tential which approaches zero for small distances. It also 
misses the proper sequence of the antisymmetric chan- 
nel to be above the symmetric one for small distances. 
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This is the result of leaving out the normalization de- 
nominator in Eq. (6). Thus these problems are corrected 
in the Heitler-London treatment. Here also the antisym- 
metric channel shows the expected Coulomb singularity 
of the hole- hole potential for small distances. In the full 
solution the energy is lowered compared to the Heitler- 
London approximation by a mutual deformation of the 
excitonic orbitals. Please note that the van der Waals 
effect is included due to the nonperturbative nature of 
the calculation. 



full calculation 
Heitler-London 
Hartree-Fock 




symmetric 



0.0 
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Figure 1: Effective interaction potentials for the strictly 
two-dimensional system: full calculation results (solid lines), 
Heitler-London approximation (dashed lines), and Hartree- 
Fock approximation (dot-dashed lines) plotted vs distance in 
units of the three-dimensional (3D) (bulk) exciton Bohr ra- 
dius a b ■ The vertical axis is given in units of the bulk exciton 
Rydberg energy Ry* . Antisymmetric channels are displayed 
in gray, while symmetric channels have black lines. 



While the considerations above are correct also in three 
dimensions, we turn now to the two-dimensional case of 
quantum structures. All vectors are to be understood 
as lying in the x-y plane, the z-axis being the growth 
direction. 

We note that all interaction potentials in Fig. 1 ap- 
proach zero from positive values and hence show a re- 
pulsion for larger distances, although this is hardly seen 
in Fig. 1. However, this behavior follows already from a 
multipole expansion of the direct potential J7d(R), which 
dominates the interaction at large distances. In such an 
approach we treat the exciton as a static charge distri- 
bution p(r, z) with cylindrical symmetry and centered at 
zero cm. coordinate. The specific form of p(r, z) depends 
on the system under investigation and will be specified 
later. Due to the charge neutrality of the exciton the 
multipole expansion starts with the dipole term oc 1/i? 3 , 
where R denotes the in plane cm. distance of the two 
excitons. Up to the quadrupole term we obtain for the 



asymptotics 



i? -> oo : U d (R) = (zf 



+ 



^Q(x 2 -z 2 ) 2 + 3(z)(z(3 2 ; 2 -z 2 ))), (10) 



where the angular brackets denote averaging over p(r, z). 
For a bulk system with spherical symmetry we see im- 
mediately that the multipole expansion vanishes as ex- 
pected. However, for the reduced symmetry of quasi- 
two-dimensional systems, we get finite multipoles also 
for in-plane circular symmetry. Even for the strictly two- 
dimensional system, where z = 0, there is a contribution 
from the quadrupole-quadrupole interaction (oc 1/i? 5 ). 
This holds as well for symmetric single QWs where 
(2) = due to the mirror symmetry along the z-axis. 
The direct interaction thus falls off as a power law rather 
than exponentially, which has been overlooked in a re- 
cent investigation of the asymptotic XX potential. 25 In 
Fig. 2 we compare the asymptotic behavior of the nu- 
merically calculated direct potential [Eq. ( A4)] , and the 
result from the multipole expansion [Eq. (10)] using the 
strictly two-dimensional charge distribution for an im- 
mobile hole p(r,z) = 5(z) [S(r) — <p 2 {r)] and the usual 
Coulomb potentials [Eq. (2)]. The multipole expansion 
holds for R > 9 0,3- The full solution shows the van der 
Waals effect and lies below the direct potential. How- 
ever, this effect is not able to overcome the repulsive na- 
ture of the asymptotic potential. To grasp this effect a 
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Figure 2: Asymptotic behavior of the potentials on a dou- 
ble logarithmic plot for the strictly two-dimensional system: 
full calculation (solid), Heitler-London (dashed) and multi- 
pole expansion (dot-dashed); symmetric states lie on top of 
antisymmetric ones. 

bit more quantitatively, we plot in Fig. 3 the difference 
between the direct potential from the approximations of 
rigid orbitals and a corresponding quantity for the full 
solution (t/^ u11 = l/2(U a + U s )). For very large distances 
i? > 10 as, we fit this difference to a van der Waals poten- 
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tial, i.e., a oc 1/R 6 power law (dashed line in Fig. 3). The 
van der Waals law holds only for distances R > 9 as . For 



10" 



■ difference: direct potential - full solution 
■van der Waals fit 
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Figure 3: The van der Waals effect in the strictly two- 
dimensional system. Deviation of the full calculation from 
the direct potential for asymptotic large distances (solid line) 
and fitted van der Waals like potential oc 1/R 6 (dashed line). 

smaller distances higher-order effects come in and spoil 
the ex 1/R 6 dependence. 

With the derived effective interaction potentials, we 
construct the two-dimensional Schrodinger equation for 
two excitons with mutual distance R in cm. space, where 
we introduce now the finite exciton mass M = m e + mh 
in the kinetic term: 



T_r_ 

M 



--R— 

RdR dR 



m 

R? 



i/*(jz)jc(«) 

= E^l(R), (11) 



where m denotes the angular quantum number. Please 
note the missing factor of 2 in the denominator of the 
kinetic term since we have to consider the reduced mass 
of two excitons, fix = M/2. 



IV. MODEL SYSTEMS 

To test the reliability of our Born-Oppenheimer-type 
method, we calculate the biexciton binding energy Bxx 
from the lowest state of Eq. (11) with m = for different 
mass ratios a — m e /mh in the strictly two-dimensional 
limit. These results are compared in Fig. 4 with vari- 
ational calculations from the literature 32 which are nu- 
merically exact. As expected, our method produces exact 
results for a = 0. For nonzero mass ratio we underesti- 
mate the binding energy slightly, e.g., for a — 0.3 we 
have an error of w 8%. Even for a = 1 the error is only 
about 12%. We conclude that the method is a reasonable 
approximation for GaAs quantum structures (a = 0.3), 
which will be under investigation in Sec. V. We plot the 



literature 

full calculation 

Heitler-London 

Hartree-Fock 




0.4 0.6 
mass ratio a 

Figure 4: Haynes factor fn = Bxx/Bx for the different po- 
tentials in the strictly two-dimensional system: full calcula- 
tion (solid line), Heitler-London potential (dashed line), and 
Hartree-Fock potential (dot-dashed line) plotted vs mass ra- 
tio a = m e /mh' The dots are numerically exact values from 
Ref.32. 



Haynes factor fjj — Bxx/Bx for the other approxima- 
tions as well (for the strictly two-dimensional system, 
Bx = ARy*). The Heitler-London approximation un- 
derestimates the biexciton binding energy significantly. 
Hartree-Fock seems to do much better, but note that for 
small a the Haynes factor is overestimated, which is in 
contrast to the variational principle. 

A simple model for a coupled quantum well structure 
is the so- called bilayer system: 39 Electrons and holes are 
confined each in infinitely narrow planes with a separa- 
tion d between the layers. This separation of unequal 
charges leads to a reduced Coulomb interaction between 
particles in different layers but leaves the potentials be- 
tween particles of the same kind unchanged, 



eh(r) 



Vr 2 + d? ' 



»(r) 



■■ v hh (r) = ^. 

r 



(12) 



The charge distribution reads now p(r, z) = 5(z)S(r) — 
5(z — d)(f> 2 (r), and the expansion [Eq. (10)] yields a finite 
dipole-dipole interaction oc 1/i? 3 , resulting in a strong 
long-range repulsion. For the many-exciton problem we 
are in particular interested in the biexciton formation. 
Fig. 5 shows the biexciton binding energy versus charge 
separation for a mass ratio of a — 0.3 and a = 0.5. 
In both cases we find a fast reduction of Bxx with in- 
creasing separation between the layers as expected. At a 
certain critical charge separation d cr it, denoted by arrows 
in Fig. 5, the biexciton ceases to exist. The same model 
system has been investigated by Tan et al. 39 using quan- 
tum Monte Carlo (QMC) technique. From their results 
at intermediate values of d, which agree nicely with our 
own calculations, the authors suggested an exponential 
decay of Bxx for large d. The inset of Fig. 5 shows the 
potential in the symmetric channel for different values of 
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the charge separation. It is clearly seen how the poten- 
tial minimum passes through zero with increasing d and 
vanishes completely if d is enlarged further (large arrow). 
Therefore a (finite) critical charge separation d cr it exists 
where the biexciton binding energy becomes zero due to 
the dominant dipole-dipole repulsion. More recently this 
critical behavior was also observed in QMC calculations 
of Lee and Needs. 40 In Fig. 6 we plot the critical charge 
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Figure 5: Biexciton binding energy in the bilayer system plot- 
ted vs charge separation at a = 0.5 (solid line) and a = 0.3 
(dashed line). The arrows denote the corresponding critical 
charge separation d cr it • Curves in the inset show the symmet- 
ric potential U a (R) for d = 0.8, 0.9, 1.0, and 1.1 in units of 
the 3D exciton Bohr radius as. 

separation vs mass ratio. Below this curve, in the shaded 
area, bound states can be formed, while above no biex- 
citons exist. For small o we obtain a rapid decrease of 
the critical charge separation. For larger mass ratios, on 
the other hand, the separation does not depend much 
on a. This behavior resembles Fig. 4, revealing a direct 
connection between biexciton binding energy and critical 
charge separation. 



V. RESULTS FOR COUPLED QUANTUM 
WELLS 

We turn now to realistic coupled quantum well struc- 
tures, having GaAs as well material and barriers made of 
Al 2; Gai_ 3 ;As. In this case we treat the Coulomb interac- 
tion v a i,(r) in single sublevel approximation, 33,41 which 
is given by 



v a b(r) = eg / dz'dz 



u 2 (z')ul(z'-z) 
yfr 2 + (z - z') 2 ' 



(13) 



Here, r is again an in-plane vector, and u e (z) and Uh(z) 
denote the confinement functions of the lowest sublevel 
for electron and hole. They enter as well the static 
charge distribution of the exciton: p(r, z) = u 2 (z)5(r) — 



no biexciton 
formation possible 




Figure 6: Critical charge separation in units of the bulk exci- 
ton Bohr radius as plotted vs mass ratio a = rn e /irih- 



u 2 (z)(j> 2 (r). Fig. 7 shows the effective interaction poten- 
tials for a GaAs/Alo.3Garj.7As coupled quantum well ge- 
ometry used by Butov 28 with a nominal (i.e., center dis- 
tance between the wells) charge separation of d = 12 nm 
(sample A in Table III). We see again the already dis- 
cussed features of the three approximation levels. The 
energy gain when going from Heitler-London to the full 
calculation yields a minimum for the symmetric chan- 
nel, which however is so weak that no biexcitons can be 
formed. This feature strongly depends on the geometry 
of the quantum wells, which is illustrated in Fig. 8. Here 
the full calculation potentials for sample A are compared 
to another one used by Snoke et al. 42 where d = 14 nm 
(sample B in Table III). Due to the larger charge separa- 
tion no minimum can be seen in the symmetric channel. 
This observation is consistent with the results obtained 
for the idealized (bilayer) model discussed in Sec. IV. 



Table III: Details of the used CQW geometry together with 
the calculated binding energy Bx of the indirect exciton. 

Sample A Sample B 



Well width L z (nm) 


8.0 


10.0 


Barrier width (nm) 


4.0 


4.0 


d (nm) 


12.0 


14.0 


d cB (nm) 


10.8 


12.7 


Static field (kV/cm) 


30.0 


36.0 


B x (meV) 


4.0 


3.5 



With these results a simple explanation for the 
regular bead pattern in the luminescence ring at 
low temperatures 27-29 has to be ruled out: It was 
speculated 25 ' 26 that the van der Waals effect could over- 
come the dipole-dipole repulsion, resulting in an attrac- 
tion between spatially indirect excitons, which would lead 
to a spontaneous patterning. Our calculation shows that 
this is not the case in agreement with recent experimental 



investigations 



43,44 



In Ref. 25 the quadrupole-quadrupole 
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Figure 7: Interaction potentials for sample A: full calcula- 
tion results (solid lines), Heitler-London potentials (dashed 
lines), and Hartree-Fock potentials (dot-dashed lines) plotted 
vs distance. Antisymmetric channels have gray lines, while 
symmetric channels have black ones. 



interaction in two-dimensional systems has been over- 
looked and hence the role of the van der Waals force 
overestimated as discussed in Sec. IV. 



0.08 



0.06 - 



"o 0.04 h 
E, 

LU 



0.02 



0.00 



sample A (d=12nm) 
sample B (d=14nm) 



25 




125 



Figure 8: Comparison of the full calculation results for sample 
A (solid lines) and sample B (dashed lines). Antisymmetric 
channels have gray lines, while symmetric channels have black 
ones. 



VI. EXCITON BLUESHIFT AND 
BROADENING 

We turn now to the calculation of the interaction in- 
duced blueshift and broadening which can be observed in 
photoluminescence experiments. 42,45,46 In a point charge 
treatment of spatially indirect excitons, we have from 



Eq. (A4) 



U d (R) = v hh (H) 



(R)-2v eh (R). 



(14) 



Plugging in Eq. (12) for the bilayer system, one gets a 
dipole-dipolc repulsion of the form 



U d (R) 



2 
R 



Vi? 2 + <P 



(15) 



Assuming a homogeneous exciton density nx , this leads 
to a blueshift, 



A = / d 2 RU d (R)n x 



d nx- 



(16) 



Since this expression is consistent with the electrostat- 
ics of a plate capacitor, it is often referred to as ca- 
pacitor formula. 47 We will derive a corrected formula 
for the blueshift and the scattering-induced broaden- 
ing using the effective interaction potentials given in 
Fig. 8. The exciton self energy is calculated in a T-matrix 
approach. 34 In the low-density limit and assuming com- 
plete spin equilibrium, we write the two-body T-matrix 
equation as 



;q|2*(z)|q">=0* - 



' 2e a ' — hz 

q 



;q'|T^z)|q" 



(17) 

The T-matrix enters the quasiparticle self-energy as 
boson-direct (D) and boson-exchange (X), 

S k ( ek )=4^{(q|T^(z)|q) + <q|T^(z)|-q)} 
q 

X " fi (£ (k + 2q)2 ^)' (18) 

where hz = h 2 q 2 jM + iO is put on shell and n.e(e) is 
the exciton distribution function, which is later taken 
as the low-density Maxwell-Boltzmann expression. The 
spin structure of the Hamiltonian yields the following 
decomposition in the limit of immobile holes: 

T D = 3T a + T s and T x = 3/2 T a - 1/2 T s . (19) 

The on-shell T-matrix needed in Eq. (18) depends exclu- 
sively on the phase shifts (q) via 48 



;q|T»(z)|±qH-5>ir- 



cot( 



»(«)) 



(20) 



We have extracted the phase shifts from the asymp- 
totics of the solution of the radial Schrodinger equation 
[Eq. (11)] for E = H 2 q 2 /M. Results for the total complex 
scattering amplitude [curly bracket in Eq. (18)] are shown 
in Fig. 9. Please note that for large momenta q, the real 
part of the scattering amplitude approaches the predic- 
tion of the capacitor formula but with an effective charge 
separation d Q g which is somewhat below the nominal one 
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Figure 9: Real (solid line) and imaginary (dashed line) parts 
of the scattering amplitude as a function of momentum for 
sample A. Included is the Maxwell distribution at T = 100 K 
to show the relevant values for q. The ordinate has been 
normalized to the capacitor value [Eq. (16)]. 



(third and forth row of Table III) . This reduction is due 
to the spatial extension of electron and hole charges along 
z. A simple argument assuming confinement wave func- 
tions for infinite barriers leads to d c s ~ d — Q.12Q1L Z . 
This is quite close to the numerical value which follows 
from replacing the point charge potential [Eq. (15)] in 
Eq. (16) by the numerically derived one. 

Using the scattering amplitude, we calculate the quasi- 
particle shift and broadening at the dispersion edge (k = 
0) and introduce correction factors ft(T), fz(T) to the 
capacitor formula 

S (0) = d— n x (h(T) - if 2 (T)) . (21) 
eoes 

The real part of this quantity is the blueshift A of the 
exciton due to the repulsive interaction, while the imag- 
inary part can be associated with a finite broadening. 
Please note that with the sign convention used, ImS is 
negative. The correction factors shown in Fig. 10 reduce 
the capacitor result dramatically. Therefore, the density 
for a measured blueshift would be underestimated by a 
factor of 10 at low temperatures. The broadening is of 
the same order of magnitude, which is consistent with 
experimental findings. 



strong depletion of the exciton gas around a given exciton 

due to the repulsive interaction. To grasp this repulsive 

correlation more directly we calculate the exciton pair- 
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Figure 10: Correction factors to the capacitor formula in de- 
pendence on temperature for sample A (solid lines) and sam- 
ple B (dashed lines). The upper panel shows f\ (blueshift), 
while the lower one shows fa (broadening). 



correlation function: 

(*t(R)^(o)^(o)* s (R) 



*J(R)*,(R))(*t,(0)*y(0) 



n s n s 



kk'q 



(22) 



VII. EXCITON-EXCITON CORRELATION 
FUNCTION 

The significant reduction of the quasiparticle shift com- 
pared to the capacitor value can be explained with a 



It has the same spin structure as the T-matrix. In the 
spin equilibrated situation investigated here, the exciton 
density n s does not depend on spin. Summing over both 
spin indices, we obtain with a partial wave decomposi- 
tion, 



8 



Ek,m eX P 


( h 2 k 2 \ 

^ Mk B T ) 


{3 \r m (R)\ 2 + \r m {R)\ 2 + (-ir (I \r m {R)\ 2 - \ \r m m 2 ) 


} 


4E k exp 


( n 2 k 2 ^ 

^ Mk B T J 





(23) 



The results shown in Fig. 11 reflect the strong repulsion 
of the excitons independent of the spin channel. It is 
interesting to compare with the pair-correlation function 
of ideal bosons and fermions having four spin degrees of 
freedom as well. Obviously in the present case, the repul- 
sive interaction between excitons is much more important 
than the bosonic nature of their statistics. 
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Figure 11: Exciton pair-correlation function vs distance at a 
temperature of T — 6K. The solid line refers to the calcu- 
lated exciton-exciton potential for sample A. The dashed line 
represents ideal bosons, while the dot-dashed line holds for 
ideal fermions. 



much richer spin structure with spatially dependent in- 
teraction potentials than with just contact interactions. 
This is due to non trivial exchange processes of the 
fcrmionic constituents. For the spatially dependence of 
the potentials, we compared three different levels of ap- 
proximation: The Hartree-Fock and the Heitler-London 
approximations as well as a newly introduced full nu- 
merical solution of the two exciton problem. We found 
a principal failure of the Hartree-Fock treatment, which 
is cured in the Heitler-London approach. The quality 
of the latter, however, turns out to be quite poor, com- 
pared to numerically exact results. With our calculated 
potentials we have investigated bilayer systems and two 
different CQWs. The charge separation d plays a fun- 
damental role: By tuning d a transition happens from 
systems with possible biexciton formation to those where 
biexcitons are not bound due to the stronger XX repul- 
sion. For two realistic CQW structures we have calcu- 
lated the quasiparticlc shift and broadening at the band 
edge which govern roughly the photoluminescence line 
shape. At low temperatures, we found a dramatic reduc- 
tion of the blueshift compared to a naive treatment of 
the CQW as plate capacitor. The broadening turns out 
to be of the same order as the blueshift. 
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Appendix A: 



In the limit of infinitely heavy holes, the interaction part of the many-exciton Hamiltonian (Table II) can be written 



as 33 



H'xx = \j dRdR';7 a (R-R')^*i(R)*l,(R')^(R')*,(R) 

ss' 

+ i J dRdR' (U a (K - R) - II s (R - R)) e ' s ') 

SS f 

*^(R)*t(R)$ s ,(R')$ s (R) _ *t( R )^t ;(R /^ s , (R ^ s(R) 

+*t,(R)*l a ,(R')*-a(R')*.(R) (1 - 25„')] • (Al) 

For the simplifications of contact potentials U P (R — R') = U P S(R — R'), our result agrees with the one derived by 
de-Leon and Laikhtman. 8 This can be seen as follows: Plugging the contact potentials into Eq. (Al) the integration 
over R' can be carried out immediately, 

H'Tx = \u a f <m£*KR)*l,(R)* S '(R)* s (R) 

+ -{U a -U s ) /'rfR^e(s-s / )*i,(R)* t _ s ,(R)*- s (R)* s (R)(l-2^). (A2) 

^ ss' 

In this limit, the third row of Eq. (Al) gives no contribution. We can also drop the Heaviside step function by inserting 
a factor of 1/2 to account for the double counting in the sum over s and s' and get 

h'xx = \u a f dR5>J(R)*^(R)*.'(R)*.(B.) 

J ss' 

+ l(U a -U') /"dR£*t,(R)*t a ,(R)*_ a (R)* s (R)(l-2<J M ,)- (A3) 

ss' 

In the language of the Hartree-Fock approximation, we can trade the symmetric and the antisymmetric potentials for 
the direct and the exchange ones. For arbitrary hole masses, Ud(R) and U X (R) have been derived in the literature. 6 
For immobile holes they reduce to the direct potential, 

U d (R) = v hh (R) + J dvdv'<f(v)v ee {r - r> 2 (r' - R) - 2 J dvv eh {v)4> 2 (v - R), (A4) 

and the exchange potential, 

U X (R) = -0 2 (R)v hh (R) - J drdrV(r)0(r - H)v ee (r - r / )0(r / )0(r / - R) + 20(R) J dr« efc (r)0(r)0(r - R). (A5) 
Taken in the contact limit U a = Ud + U x and II s = Ud — U x , Eq. (A3) can be written as 

H'x X = \vdl dRX;*J(R)*J'(R)*.'(R)*-(R) 

ss' 

+- A u x ydR^[2*t( R )^t /(R) ^ s , (R) ^ s(R) + ^t /(R) ^t_ s/(R) ^_ s(R) ^ (R)(1 _ 25ss , ) ] _ (A6) 

ss' 

The Fourier transform of this result yields the Hamiltonian derived in Ref. 8. 
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